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Abstract 

We report the results of mean field and the Monte Carlo study 
of the dynamic magnetization-reversal transition in the Ising model, 
brought about by the application of an external field pulse applied 
in opposition to the existing order before the application of the pulse. 
The transition occurs at a temperature T below the static critical tem- 
perature Tc without any external field. The transition occurs when 
the system, perturbed by the external field pulse competing with the 
existing order, jumps from one minimum of free energy to the other 
after the withdrawal of the pulse. The parameters controlling the 
transition are the strength hp and the duration At of the pulse. In 
the mean field case, approximate analytical expression is obtained for 
the phase boundary which agrees well with that obtained numerically 
in the small At and large T limit. The order parameter of the tran- 
sition has been identified and is observed to vary continuously near 
the transition. The order parameter exponent P was estimated both 
for the mean field (/? = 1) and the Monte Carlo (/3 = 0.90 ± 0.02 in 
two dimension) cases. The transition shows a "critical slowing-down" 
type behaviour near the phase boundary with diverging relaxation 
time. The divergence was found to be logarithmic in the mean field 
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case and exponential in the Monte Carlo case. The finite size scaling 
technique was employed to estimate the correlation length exponent 
V (= 1.5 ± 0.3 in two dimension) in the Monte Carlo case. 

1 Introduction 

The dynamic response of pure Ising systems to time dependent magnetic 
fields is being studied intensively these days (see e.g. |]l], H and references 
therein). In particular, the response of Ising systems to pulsed fields have 
recently been investigated §]. The pulse can be either "positive" or "neg- 
ative". At temperatures T below the critical temperature Tc of the corre- 
sponding static case (without any external field), majority of the spins orient 
themselves to a particular direction giving rise to the prevalent order. If the 
external field pulse is applied along the direction of the existing order, it is 
called a positive pulse and if the pulse is applied opposite to the existing 
order, it is called a negative pulse. The effect of positive pulse has been 
studied by Acharyya et al [Q, whereas the occurrence of a magnetization- 
reversal transition as a result of the application of the negative pulse has 
been reported in an earlier work 0]. We report here the results of detailed 
investigation of this dynamic magnetization-reversal transition for pure Ising 
models under a pulsed field. 

In the absence of any symmetry breaking field, for temperatures below the 
critical temperature of the corresponding static case (T < Tc), there are two 
equivalent free energy minima with average magnetizations +mo and —mo. If 
in the ordered state the equilibrium magnetization is +mo (say) and the pulse 
is applied in the direction opposite to the existing order, then temporarily 
during the pulse period the free energy minimum with magnetization —vhq 
will be brought down compared to that with +mo. If this asymmetry is 
made permanent, then any nonzero field (strength), which is responsible for 
the asymmetry, would eventually induce a transition from +mo to —mo. 
Instead, if the field is applied in the form of a pulse, the asymmetry in the 
free energy wells is removed after a finite period of time. In that case, the 
point of interest lies in the combination of the pulse height or strength [h^ 
and its width or duration (At) that can give rise to the transition from 
-|-mo to — mg. We call this a magnetization- reversal transition. A crucial 
point about the transition is that it is not necessary that the system attains 
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its final equilibrium magnetization —mo during the presence of the pulse; 
the combination of hp and At should be such that the final equilibrium 
state is attained at any subsequent time, even long time after the pulse is 
withdrawn. The phase boundary, giving the minimal combination of hp and 
At necessary for the transition, depends on the temperature. As T — Tc, 
the magnetization reversal transition occurs at lower values of hp and/or At 
and the transition disappears at T > T^. 

In this paper we have given both the mean field (MF) and the Monte Carlo 
(MC) results for the transition. The MC studies have been carried out for a 
two dimensional {d = 2) lattice of Ising spins. The phase boundaries in the 
hp — At plane (each for a fixed T) are obtained for both the MF and the MC 
cases. Approximate analytical expressions are also obtained and compared to 
these phase boundaries. The order parameter (O) for the dynamic transition 
has been identified and at a fixed T its variation with the driving parameters 
hp and At has been studied. The observed continuous variation of O indicates 
the nature of the transition to be comparable to the second order type static 
transitions. The critical exponents for the order parameter variations near 
the phase boundary has been obtained for both the MF and MC cases. We 
also employed the finite size scaling method (see e.g 0) to estimate the 
correlation length exponent u for the transition in the MC studies. We 
observe significant "critical slowing down" near the phase boundary and the 
behaviour of the relaxation time r has been studied in both the cases. In 
the MF case, we have obtained an approximate analytical expression for r 
indicating clearly a different kind of divergence of r as compared to that in 
the static case. 



2 Model 

We have taken Ising Model for both the numerical simulation and mean 
field study. The Hamiltonian of nearest neighbour Ising system without any 
disorder is 

h = -jY1 s,Sj-h{t)J2Si, (1) 

<ij> i 

where Si = ±1 represents the Ising spins at lattice site i and J denotes 
the nearest neighbour interaction strength. The time dependent external 



3 



magnetic field h{t) is applied in the form of a pulse of duration At 

h{t) = -hp, for to <t <to + At 

— 0, otherwise. (2) 

to is taken to be much larger than the relaxation time of the unperturbed 
system so that the system is guaranteed to reach a state of equilibrium before 
the pulse is applied. The average magnetization m is given by < Si >, where 
the angular brackets represent thermal average. By the time t = to, the 
system has reached its equilibrium state with the magnetization m(t) = ±mo. 
This is the state before the application of the pulse, which is applied at 
time t — to to compete with the initial magnetization. We then want to 
look at the dynamics of the system under a pulsed field starting with the 
initial condition m(to) = (say). Depending on the strength (hp) or 

the duration (At) of the pulse, the system has two choices after the pulse 
is withdrawn: it can either go back to the original ordered state {m{t — 
oo) = +mo), or it can switch to the other equivalent equilibrium ordered 
state {m{t = oo) = —mo). Fig. 1 shows schematically these behaviours. 
The result, at any finite temperature below T^, naturally depends on the 
strength and the duration of the pulse. SpecificaUy, we observe that the 
result depends on rriyj = m{to + At), the average magnetization at the time 
of withdrawal of the field. The sign of m^, which in turn depends on the 
combination of hp and At, governs the transition. | \ was found out to be 
the appropriate candidate for the order parameter: if rriy, becomes negative, 
on an average one observes a magnetization-reversal transition; however if m^^ 
is positive, the system is most likely to return back to its original equilibrium 
state. In fact, it is guaranteed to be so in the mean field case. However in 
MC simulations, there exists occasional fluctuations due to which the system 
may finally arrive at the +mo state starting with negative and vice versa. 
We have identified the hp — At phase boundary at a particular temperature 
which gives the optimum combination of the driving parameters {hp and At) 
that can force the system to a final state with magnetization —mo, starting 
from an initial state with magnetization +mo or vice versa. We define the 
relaxation time r as the time taken by the system to reach its final equilibrium 
state from the time of withdrawal of the field (at t = to + At). The relaxation 
time increases as the value of | m^ | approaches zero or equivalently as one 
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approaches the phase boundary from either side. In Fig. 2 the typical MC 
results (on a square lattice) show how the transition can be brought about by 
either increasing hp (see Fig. 2(a)) or At (see Fig. 2(b)). One can also note 
from these figures how the relaxation time r increases as one approaches the 
phase boundary. 



3 Mean Field Study 

The mean field equation of motion for the average magnetization m{t) of the 
system is 

dm , f m + hit) \ 

_._,„ + ta„h^^iij. (3) 

where h{t)is given by (|^). Here we have assumed Jn = 1, where n is the lat- 
tice coordination number. With this choice, the critical temperature in the 
static limit (T^) becomes unity. All the mean field calculations are performed, 
therefore, at a temperature T < 1. The equation was solved numerically to 
obtain the phase boundaries, shown in Fig. 3. A point on a particular 
phase boundary gives the optimal combination of hp and At that can in- 
duce the transition from a state with magnetization +mo to —mo, where 
mo = tanh(mo/T). The axes side of the boundaries correspond to the return 
to original equilibrium state, whereas one gets a magnetization-reversal tran- 
sition for combinations of hp and At beyond the phase boundary. Because 
of the absence of any fluctuations in the mean field equation (^, there exists 
a finite coercive field for T < and therefore one cannot bring about the 
transition just by increasing At if hp does not exceed the coercive field value. 
Hence the phase boundary becomes parallel to At axis for large values of At. 

For large T, the mean field phase boundary can be estimated approxi- 
mately by solving the linearized mean field equation 

dm hit) 1 — T , . 

_ = -em + ; e = — — , (4) 
dt T T ^ ^ 

where hit) is given by (Q). The solution of (|^) for t > to can be written as 
m(t) = Uo - ^ J exp[e(t - to)] + ^. (5) 
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It may be noted from the solution that t has to be close to to in order 
to keep the value of m{t) small, so that the linearization in is valid. 
The magnetization-reversal transition occurs if m(to + At) < 0. The phase 
boundary can therefore be obtained from (|) by putting = m(to + At) = 
0. The resulting equation for the phase boundary is then 

P l-exp(-eAt)- 

In Fig. 3, this analytic result for the phase boundary has been compared 
with those obtained by solving numerically equation (|^). As one can clearly 
see from the figure that the agreement is good in the small At region of the 
phase boundaries for large values of T. 

As mentioned before, the magnetization m^; at the time of withdrawal 
of the pulse, seems to be the crucial quantity governing the transition. The 
sign of solely decides the final equilibrium state of the system out of the 
two equilibrium choices. Therefore we define the mean field order parameter 
(O) as following 

O = m^6{mui), (7) 
where the step function 6 is defined as 



6{x) = 1, for X > 
= 0, otherwise. 

The nature of variation of O with | /ip — /i^ | at different values of At and T is 
shown in Fig. 4. Here = /ip(At,T) is obtained from the phase boundary. 
We fitted the order parameter variations to the power law form 

o ~| /ip - h; 1^ (8) 

and found the value for the exponent /5 ~ 1 by fitting the numerical results 
for different values of At and T. 

The relaxation time (r) grows as the phase boundary is approached from 
either side and shows a divergence on the boundary. We measure the re- 
laxation time by measuring the time required by m(t) to reach the final 
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equilibrium value ±mo, with an accuracy of O(10~^), from the time of with- 
drawal of the pulse. According to and (H), for a fixed At, ~| hp — h^ \. 
Therefore r is also expected to diverge as vanishes. Fig. 5 shows that 
this is indeed the case and the growth of the relaxation time was found to 
be logarithmic in nature. That the relaxation time will diverge logarith- 
mically at the phase boundary at any temperature below the static critical 
temperature (T < 1) can be shown analytically. Let us follow the mean field 
dynamics of the system after the withdrawal of the field. The system starts 
with a magnetization and evolves according to (|^) with h{t) = to reach 
the final state of equilibrium characterized by magnetization ±mo. Keeping 
the value of m/T small we expand tanh of (|^) upto the cubic term : 

dm o 

— — = em + am , (9j 

where a = —1/ (3T^). Thus we can write 

, r^o dm 
dt 



Jm^ em + am^ 
or, 

r = -i^ + i^^^P^ + C(T), (10) 

where C{T) is a constant depending on temperature only. Now, if the com- 
bination of hp and At is such that one starts with very small value of m^, 
then one gets 

r ~ — In I | . (11) 

It can be easily verified that even if the O(m^) term in the expansion of (|^) 
is kept, the final result (|TTD does not get modified in the m^, limit. 
The above form for the variation of r fits accurately the numerical results : 
at T = 0.6, 1/e = 1.5, and this matches exactly with the numerical result 
as shown in Fig. 5. The above analysis also makes it very clear that the 
logarithmic divergence of m^, implies a similar divergence with \hp — hp\ for 
a fixed At. It may be mentioned here that for the dynamic magnetization- 
reversal transition, the above divergence in r occurs at any T < T^.. The 
same equation ( |ri|) also gives the well known | T — divergence of r 

for the static transition at T = = 1. It may also be noted from ([11]) that 
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the logarithmic dependence of t on \ hp — \, through its dependence on 
rriu], for this dynamic transition at hp{At, T < Tc) is quahtatively different 
from the divergence of r (~| T — T^ |^^) for the static transition in the same 
model (equation |^) at T = Tc. 

4 Monte Carlo Study 

For the MC simulation we have taken Ising spins on a square lattice of size 
L X L, with L = 200 for typical studies. Each spin of the lattice were 
updated sequentially using the Glauber single spin flip dynamics [|]. Phase 
boundaries at different temperatures below the static critical temperature 
(Tc ^ 2.27) shows qualitatively similar nature to that in the mean field case. 
However, unlike the mean field phase diagrams, the presence of fluctuations 
causes the MC phase diagrams to touch the abscissa asymptotically. When 
the contributions of fluctuations become important at higher values of T and 
small values of hp, the above MF theory fails. If hp (as in the large 
At region of the phase boundary), one can use the picture of nucleation of 
a single domain. The classical nucleation theory of Becker and Doring (see 
e.g. 10]) suggests that the nucleation rate I ~ exp[— T(Zc)/T] is given by the 
optimality condition Ic = [a{d — l)/2d\/hp of free energy F{1) = 2hpl'^ + al'^~^ 
for the formation of a droplet or domain of linear size / under field hp. Here a 
is proportional to the surface tension for the formation of a droplet. Equating 
the growth rate given by Becker and Doring nucleation rate / with the inverse 
pulse width, one gets At ~ exp(l//ip~^), suggesting 

/iplnAt = constant, (12) 

along the phase boundary in two dimensions. It agrees fairly well with the 
MC estimated phase diagram in the low hp (i.e. large At) limit 0, ^. For 
large values of hp, the critical droplet size Ic ^ 1/hp being much smaller than 
the system size L, many droplets grow simultaneously and the transition 
occurs due to coalescence of the droplets and not by the above-mentioned 
process of growth of a single droplet. The equation (|T^), therefore, fails to 
fit the phase boundary for large values of hp and hence for small values of 
At. 

Similar to the MF case (|^), we define the order parameter in the same 
way : O = mw6{mw), where 6 is the step function and ruw is the average 
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magnetization at the time of withdrawal of the pulse. Fig. 6 shows the 
variation of this order parameter with | hp — hp \ for different values of At and 
T. Typical number of initial seed values taken to generate the configurations 
over which each data point is averaged is 500. The variation of O was again 
fitted to the power law form with the order parameter exponent f3 and 
we find /? ~ 1 by fitting the MC results for four different values of At and 
four different values of T. 

Here also we estimate the relaxation time r by measuring the time (MC 
steps) required for m{t) to reach the final equilibrium value ±mo from the 
time of withdrawal of the pulse with a predefined accuracy O(10~^). Again 
T was found to diverge as the phase boundary (at any fixed T) is approached 
from either side. Fig. 7 shows that the divergence of r occurs at the point 
where m^^, vanishes at the phase boundary. It is not possible to average the 
relaxation time for different realizations of the dynamics for a particular m^, 
as different realizations produce different values of rriy^. Therefore a small 
range of was taken, instead of a particular value, to take an average of 
r. A typical value of the range of rriw over which averaging has been done 
is 0.01. It was observed that r depends on the driving parameters hp, At 
and temperature T only implicitly through the quantity rriy^. The divergence 
of T was found to be very sharp and it fitted with an exponential function 
(r ~ exp I rriw |) unlike the mean-field case where r was seen to behave 
logarithmically with | m^, |. 

We fitted above the MC results for L = 200 and extracted the value of 
the exponent /3, assuming the finite size effects to be negligible. As hp —>■ 
hp, for finite L one should consider the effective hp to be a function of L. 
Therefore one should consider an additional scaling function of L/C, in (|^), 
where ^ ~| /ip — /i^ |~'' is the correlation length with exponent u. Assuming 
that ^ ~ L for such MC cases, we can write the appropriately modified form 
of (0) as 

Or^L-^/''f[{hp-h;)L'^'']. (13) 

We have fitted the data for different L to the scaling form (0). The data for 
L= 50, 100, 200 and 400 for a particular value of At ( =5) and two different 
values of T were scaled to obtain the fitting values of the exponents as shown 
in Fig. 8. Typical number of averages over different initial seeds taken for 
the MC data was 2500 for L = 50 and 50 for L = 400. With such a scaling 
fit procedure, hp could be obtained with an accuracy of 0.001. This gives 
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(3 = 0.90 ± 0.02 and z/ = 1.5±0.3. It was observed that the quahty of fitting 
does not change appreciably for variation of the fitting parameters within 
the error limits specified above. 

5 Discussions 

The recently reported 0] dynamic magnetization-reversal transition in pure 
Ising model under a pulsed (negative) magnetic field has been studied exten- 
sively in this paper employing the mean field approximation as well as the 
Monte Carlo technique in two dimension. Accurate estimates have been made 
for the phase boundary hp{At,T), both in MF and MC cases. Approximate 
analytic estimates of the phase boundary have been made in both the cases 
and compared with the numerical results obtained. The order parameter for 
this dynamic transition has been identified and it is given by the average 
magnetization at the time of withdrawal of the pulse. In the mean field ap- 
proximation the order parameter was found to vary linearly with | hp — \, 
indicating the order parameter exponent (3 = 1; while in the MC case /3 was 
found to be around 0.9 for d = 2. The application of the finite size scaling 
method gives the value of the critical exponent u for the correlation length 
to be nearly 1.5 for d = 2. Besides the length scale, a time scale was also 
found to diverge both in case of mean field and the Monte Carlo studies. In 
the mean field case, the relaxation time r was found to diverge logarithmi- 
cally with m^, or equivalently with | hp — |, as the phase boundary is 
approached. This has been demonstrated to be true analytically also. How- 
ever, the divergence of r was found to be stronger in the Monte Carlo case, 
where r diverges exponentially with | | or with | hp — hp \^ (with (3 ~ 0.90 
in d = 2) as the phase boundary is approached. The finite size scaling results 
suggested that the correlation length diverges as \ hp — hp at the phase 
boundary with i/ ~ 1.5 in d=2. The existence of the divergence of both the 
length scale and the time scale indicates the thermodynamic nature of this 
intriguing dynamic magnetization-reversal transition in the Ising model. 
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Figure Captions 

Figure 1. Schematic time variation of the external field h{t) and of the cor- 
responding response magnetization m{t). The thin line (smaller hp) shows 
the case of no-transition, whereas the thicker line (larger hp) shows the oc- 
currence of transition. The various quantities studied in the paper like the 
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relaxation time r, magnetization at the time of withdrawal of the field m^, 
etc are also indicated. 

Figure 2. Monte Carlo results for the time variation of the magnetization 
m(t) (dotted hne) against the time variation of the field h{t) (solid line) at 
T = 1.0 for (a) hp = 1.04, At = 10; (b) hp = 1.04, At = 15; (c) hp = 1.04, 
At = 28; (d) hp = 1.04, At = 10; (e) hp = 1.17, At = 10; (f) hp = 1.34, 
At = 10. 

Figure 3. Mean field phase boundary in the hp — At plane for three different 
values of T. The dotted lines correspond to the numerical solution of equation 
(^ while the solid lines give the corresponding theoretical prediction (from 
equation (|n|)). 

Figure 4. Log-log plot of the variation of order parameter O against | 
hp-h^p\ for At = 1, T = 0.9 (O); At = 1, T = 0.6 (A); At = 2, T = 0.8 
(□); At = 5, T = 0.9 (+); At = 10, T = 0.8 (x) and At = 10, T = 0.6 (*) 
in the mean field case . The h^ values are obtained from the corresponding 
phase boundaries (Fig. 3). 

Figure 5. Divergence of the relaxation time r in MP case. The solid line 
shows the logarithmic fit to the numerical solution of @. 

Figure 6. Log-log plot of the variation of order parameter O against | 
hp-h'j,] for At = 2, T = 2.0 (O); At = 1, T = 1.5 (+); At = 5, T = 1.5 
(A); At = 10, T = 2.0 (x); At = 5, T = 1.0 (□) and At = 1, T = 1.0 (^) 
in the MC case. The hp values are obtained from the corresponding phase 
boundaries in the MC case. 

Figure 7. Divergence of the relaxation time r for (a) T = 1.5; At = 3 (O), 
At = 5 (+) and (b) T = 2.0; At = 1 (O), At = 3 (+) in MC case. The solid 
lines indicates the exponential fit to the numerical data. 

Figure 8. Finite size scaling analysis of the MC data: 0{L)/L^/^ plotted 
against \hp-hl \ L^l"" for a fixed At. (a) T = 1.5; L = lOO(O), L = 200(+), 
L = 400(n); (b) T = 2.0; L = 50(0), L = 100(+), L = 200(n), L = 400(x). 
The best collapsed data are shown with fitting values /3 = 0.90, u = 1.5 and 
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1.088 in (a) and M = 0.720 in (b). 
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